Universal Quantum Localizing Transition of a Partial Barrier in a Chaotic Sea 
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Generic 2D Hamiltonian systems possess partial barriers in their chaotic phase space that restrict 
classical transport. Quantum mechanically the transport is suppressed if Planck's constant h is large 
compared to the classical flux, h ^> <£, such that wave packets and states are localized. In contrast, 
classical transport is mimicked for h <C <3>. Designing a quantum map with an isolated partial barrier 
of controllable flux $ is the key to investigating the transition from this form of quantum localization 
to mimicking classical transport. It is observed that quantum transport follows a universal transition 
curve as a function of the expected scaling parameter <$>/h. We find this curve to be symmetric to 
<&/h = 1, having a width of two orders of magnitude in and exhibiting no quantized steps. 

We establish the relevance of local coupling, improving on previous random matrix models relying 
on global coupling. It turns out that a phenomenological 2 x 2-model gives an accurate analytical 
description of the transition curve. 

PACS numbers: 05.45.Mt, 03.65.Sq 



In the phase space of generic two-degree-of-freedom 
(2D) Hamiltonian systems regions of regular and chaotic 
motion are dynamically separated by impenetrable bar- 
riers. Within a chaotic region so-called partial barriers 
are ubiquitous. They divide it into distinct sub-regions, 
connected by the turnstile mechanism, which works like a 
revolving door between two rooms. The volume in phase 
space, which is transported across the partial barrier in 
each direction per time is the flux <I>. Partial barriers can 
originate [1] from a cantorus or the combination of the 
stable and unstable manifold of a hyperbolic fixed point. 
A hierarchy of these partial barriers gives rise to a power- 
law decay of correlations and of Poincare recurrence time 
distributions [2J. 

What is the implication of a partial barrier on the cor- 
responding quantum system? In 1984 MacKay, Meiss, 
and Percival [3 conjectured that the flux <£, an area 
in phase space, has to be compared with the size h of 
a Planck cell to judge the quantum implications. For 
h ^> quantum transport is suppressed, while for h <C $ 
classical transport is mimicked [U{8]. Thus, the existence 
of a partial barrier in the corresponding classical system 
can be conceptualized as being responsible for partially 
localizing the quantum dynamics. As is well known, but 
still remarkable, quantum mechanics allows for both the 
suppression or enhancement of transport through local- 
ization [9j [10] or tunneling phenomena, respectively. 

Alternatively one can understand the suppression of 
transport in the time domain, where one has the Heisen- 
berg time th and the dwell time with th/t^ ~ &/h. 
For th < Td a typical classical orbit of length up to th 
does not cross the partial barrier and stays in the ini- 
tial region. To the extent that the basic semiclassical 
theory is valid (neglecting tunneling and diffraction, for 
example), the properties of the quantum system are de- 



termined by such orbits, and quantum transport must be 
suppressed. 

The quantum suppression of transport for h ^> has 
consequences for the time evolution of a localized wave 
packet initially associated with a phase space region on 
one side of the partial barrier. It cannot acquire a sub- 
stantial weight on the other side of the partial barrier 
even in the limit of arbitrarily large times; see Fig. [I] left 
inset. This is reflected in the eigenstates having much 
stronger projection in either one of the two sides; see 
Fig. [3] left insets. In contrast, for h <C $ wave packets in 
the long-time limit as well as eigenstates extend to both 
regions as if the partial barrier were not present. This 
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FIG. 1. (Color online) Asymptotic transmitted weight Woo 
vs. $>/h for <I> « (pluses), (crosses), 7^0 (circles), 

^ = IM'---'6loo com P are d to Eq. Q (solid line). Insets: 
Husimi representation of time-evolved wave packet for $ « 
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corresponds to the classical behavior where in the long- 
time limit chaotic orbits explore both regions ergodically. 

The quantum localizing transition between quantum 
suppression and classical transport has been studied the- 
oretically and experimentally for multiphoton ionization 
of atoms [6], cesium atoms in optical potentials [11 , and 
microcavity lasers [12]. The most extensive study goes 
back to Bohigas, Tomsovic, and Ullmo (BTU) on coupled 
quartic oscillators where seven chaotic regions are sepa- 
rated by six partial barriers [7] . They model the quantum 
mechanism of a partial barrier by globally coupled ran- 
dom matrices with a transition parameter determined by 
<&/h, which were previously used to describe symmetry 
breaking [13] . They found good agreement for the impli- 
cations of partial barriers on spectral statistics and wave 
packet dynamics. 

While the quantum localizing transition is partially un- 
derstood, the full quantitative transition even for an iso- 
lated partial barrier so far is not. In particular one is 
interested in the transition curve, including its universal 
scaling, center, width, shape, and whether it has quan- 
tized steps. This is a prerequisite for understanding the 
critical essence in quantum mechanics of the turnstile. 

In this paper we analyze the quantum localizing tran- 
sition with the help of a designed quantum map with an 
isolated partial barrier of controllable flux <£. By study- 
ing wave packet dynamics and eigenstate properties, it is 
found that the transition from quantum suppression to 
classical transport is universal with the expected scaling 
parameter is symmetric to &/h = 1, exhibits no 

quantized steps, and has a 10%-90% width of two orders 
of magnitude in <I>//i; for simplicity the results presented 
are quoted for a system with only two regions, both of 
comparable phase space areas. It is shown that the criti- 
cal essence in quantum mechanics of a turnstile is a local 
coupling mechanism in contrast to the previously used 
global coupling scheme. We give an analytical descrip- 
tion of the universal transition based on a phenomeno- 
logical 2 x 2-model. The findings are confirmed with the 
more familiar standard map, which plays an important 
role in studies of quantum chaos and localization [TO] . 

Consider a family of designed area preserving maps 
F of the two-torus, see Fig. [2] for an illustration. The 
phase space consists of a large chaotic sea between 
the regular tori. There is a hyperbolic fixed point at 
(QiP) = (|jPfix)j whose stable and unstable manifold 
can be used to construct a partial barrier separating 
the regions 1 and 2 of approximately the same size, 
A\ w A2. The region between the partial barrier and 
its preimage defines the turnstile areas of size pQ. 
The map F was designed such that this partial barrier 
is well isolated with a small tunable flux. This is 
achieved by a composition of two maps, F = F rot o Fki c k- 
Here F^ c ^ originates from a kicked Hamiltonian and 
is given by (q',p f ) = (q + T'{p%p* - V f (q f )/2) using 
p* = p - V'(q)/2 with V'(q) = sin(27r<2) /(4tt) and 



(a) (b) 




FIG. 2. (Color online) (a) T'(p) of the designed map F with 
$ ~ 2^0 ' 3 ) Corresponding phase space with a large chaotic 
sea (dots) between the regular tori (lines) and zoom with hy- 
perbolic fixed point (black circle), partial barrier (solid green 
line) separating regions 1 and 2, its preimage (dotted green 
line) and the turnstile areas of size 



piecewise linear T'(p), see Fig. [2|a). In order to destroy 
additional partial barriers not related to the fixed 
point we use a map F rot , which rotates points inside a 
circular region in phase space specified by center, radius 
and rotation angle while points outside this region are 
unchanged [18]. The parameters of Fkick and F rot for 
the three considered cases of F with fluxes $ ~ 3^0? 
g^, and ^0 are gi yen i n Refs. [14] [15]. Quantum me- 
chanically the system is described by a unitary operator 
U = ^ r ot^kick acting on a Hilbert space of finite size TV 
with effective Planck's constant h = 1/N. Here Z7ki c k = 
exp{-in<7)/(2fi)} exp{-iT(p)/ft} exp{-i^)/(2ft)} 
and U rot = U ho P ho + (1 - P ho ) where P ho is a projec- 
tor built from harmonic oscillator eigenstates within 
the rotating region and U^o gives the time evolution 
corresponding to the classical rotation angle [T8] . 
The time evolution of a wave packet ip(t) is given by 
i/;(t+i) = Uil>(t). 

In order to quantify the quantum transition of a partial 
barrier we define the (relative) asymptotic transmitted 
weight of a wave packet ipi started in region 1 as 

w<x#i] : = d • I 1 ) 

It is the time-averaged transmitted weight divided by 
the corresponding classical weight /i^ . As initial state 
we choose po = 0.7 and as a measure hb 2 the probability 
of ^i(t) within the region B 2 (p G [0.175,0.325]). The 
classical weight /i^ 2 is the relative time a long chaotic 
orbit spends in Assuming ergodicity in the chaotic 
sea of size it can be determined by /i^ = B2/A c \ l . 
The definition of Woo implies that it makes a transition 
from for h ^> & to 1 for h <C In the latter case 
this happens because any initial state becomes uniformly 
distributed at large times as would a classical distribution 
of trajectories. 

Figure [l] shows the resulting vs. <&/h for three 
different fluxes averaged over 100 values of the Bloch 
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phase and 100 time steps after time T = 2 20 . All data 
sets fall on top of each other under this scaling, i.e. <f>/h is 
indeed the correct scaling parameter. We expect that the 
quantum localizing transition for any partial barrier fol- 
lows the same universal curve as a function of $>/h. This 
expectation assumes that the relative volumes of regions 
1 and 2 are of the same order and that the mixing time 
is much shorter than the dwell time. Figure [T] shows that 
on a logarithmic scale the transition is symmetric with 
respect to the point $>/h = 1, Woo = |. Thus the transi- 
tion point is reached when flux <I> and Planck's constant 
h are equal. The transition is found to be smooth with no 
indications for quantized steps at integer values of $/h. 
By defining the transition region as the interval in <$>/h 
for which Woo G [0.1,0.9], the width is seen to be two 
orders of magnitude, indicated by the arrow in Fig. [T] 
The overall behavior of the transition is well described 
by the symmetric curve of Eq. Q, which results from 
a phenomenological 2x2 matrix model, see below. For 
the smallest &/h (below 0.05) there is an upward trend 
of Woo compared to the symmetric curve. We attribute 
this deviation to the effect of tunneling across the en- 
tire partial barrier, which occurs in addition to quantum 
transport through the turnstile region. This is left for 
future investigation. 

Complementary to the time evolution of wave packets, 
consider properties of eigenstates of the quantum map. 
States are either contained in region 1 or in region 2 in 
the case of h ^> <£, see Fig. [3] left inset. For h <C $ the 
states extend over the whole chaotic region, see Fig.[3]top 
right inset. In addition there exist regular eigenstates 
localized on the invariant tori and scarred states, e.g. 
on the hyperbolic fixed point. We define the average 
eigenstate equipartition measure (^12) by 
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as the sum of the products of the relative weights of eigen- 
states in each measuring box B\ (p G [0.675,0.825]) 
and B2 (p G [0.175,0.325]) contained in regions 1 and 2, 
resp. The relative weight compares the measure /jLbA^j] 
with the corresponding classical weight /i^.. Eigenstates 
which are almost zero in one of the regions give no con- 
tribution, while the contribution is 1 for eigenstates, 
which are almost equipartitioned. The prefactor l/iV c h 
is chosen such that in the semiclassical limit (^12) ap- 
proaches 1, as in this limit all iV c h = A^/h chaotic states 
contribute 1, while regular states give no contribution. 
Figure [3] shows (1V12) for three different fluxes <£. One 
observes the same transitional behavior as for Woo , again 
well described by Eq. Q. In fact, one can show that 
(11)12) = (woolipi]}^ if Woo is averaged over initial states 
ipi which form an orthonormal basis in B\. 

In order to phenomenologically describe the transition 



we propose a unitary 2x2 matrix model 
U 



y/l — V 2 V 



-y/l — v : 



(3) 



Here the deterministic variable v G [0, 1] describes the 
turnstile coupling between two sites representing the 
chaotic regions A\ and A^. Following from unitarity the 
diagonal entries have magnitude \J\ — v 2 . The lower en- 
try has a minus sign such that for v = the eigenval- 
ues are not degenerate. The eigenvalues are ±1 inde- 
pendent of v and the normalized eigenvectors are <j)± = 
l/v // 2(±c±, v/c±) with c± = y/l ± y/l — v 2 . According 
to Eq. ([2| the average eigenstate equipartition measure 
is (^12) = v 2 , where the quantum measures n B .[4>±\ are 
given by the squared z-th element of the eigenvectors </>± 
and the classical expectations are fig = /i^ = |. For 
the asymptotic transmitted weight Woo, Eq. ([!]), i.e. for 
a wave packet started on one site and measured on the 
other site, we find the same result, Woo = v 2 . The pa- 
rameter v of this model can be related to the scaling pa- 
rameter &/h of quantum maps with a partial barrier: We 
identify the transmission probability v 2 with the relative 
classical flux ^/A^, i.e. v 2 = 3>/Aj. Moreover, we choose 
the Planck cell associated with each site of the 2x2- 
model to be the sub-region of Ai that is not transmitted, 
h = Ai — This choice makes the regions associated 
with h and <£> disjoint and thereby allows for arbitrary 
ratios of <fr/h. This finally gives for the 2 x 2-matrix 
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FIG. 3. (Color online) Average eigenstate equipartition mea- 
sure (W12) for the designed map (gray crosses, same parame- 
ters as in Fig. [I]) compared to Eq. Q (solid line), the BTU 
model (gray dashed line) and the channel coupling model with 
discrete (orange circles) and continuous transmissions (light 
orange line) using 1000 realizations with Ni — N2 = 500. In- 
set: Husimi representation of typical eigenstates for $ w -55, 
h = (left) and h — (top right). Illustration of the 

matrix structure of Uo and U c , where black lines represent 
unit entries (bottom right). 



4 



model 

-oo = (w 12 ) = ITWX . (4) 

Quite amazingly this phenomenological model gives a 
very good description of the transitional behavior of the 
map data for the entire range from $>/h <C 1 to $>/h ^> 1 
(apart from the deviation for $/h < 0.05 attributed to 
tunneling at small &/ti), see Figs. [I] and [3] This is remi- 
niscent of the success of the Wigner surmise using 2x2 
matrices to describe universal spectral statistics. As no 
system specific properties were used in the derivation 
of Eq. Q, except for the scaling parameter 3>/ft, this 
gives further support for the universality of the transi- 
tion curve. We expect that the universality also extends 
to time continuous Hamiltonian systems like billiards. 

In order to get an insight into the quantum mech- 
anism of a partial barrier we now study appropriately 
adapted random matrix models. In the BTU model [7] 
two matrices of the Gaussian orthogonal ensemble, rep- 
resenting two chaotic regions, are globally coupled with 
a strength determined by $>/h. Figure [3] shows that for 
$>/h < 10 this model overestimates the value of (W12) 
found for the quantum map. We attribute this discrep- 
ancy to the global coupling of the BTU model. Instead 
we propose to model the classical turnstile mechanism 
by a local coupling via a channel. In a unitary model 
U = Uq • U c we decompose the dynamics into a coupling 
matrix U c modeling the turnstile transport multiplied by 
an uncoupled matrix Uq modeling the mixing in each of 
the regions 1 and 2, see Fig. [3] inset. The coupling matrix 
U c is an identity matrix, where the central 2n x 2n block 
has ones on the anti-diagonal. It couples the two regions 
via n modes for each direction of the channel. This mod- 
els the directed transport of a turnstile in the classical 
system. The matrix Uq is block diagonal consisting of 
two matrices of the circular orthogonal ensemble of sizes 
N/2 x N/2 and the limit of large N is considered. This 
model has only one parameter, namely the number of 
modes n = $>/h. The resulting (W12) for this unitary 
channel coupling model is shown in Fig. [3] and is in very 
good agreement with the numerical data of the map F. 
It is stressed that this agreement has been obtained with- 
out any fitting parameter. The model can be extended 
to non-integer values of &/h by the continuous transmis- 
sions of modes propagating through a channel [14 , see 
Fig. 3. We observe that these models with local cou- 
pling are better describing the data for the map F for 
&/h < 10 compared to the global coupling BTU model 
and on the same level as the phenomenological 2x2- 
model. This suggests that the turnstile mechanism of 
the classical system is quantum mechanically described 
by a local coupling via a channel. 

We now show results for the generic standard map, 
(q',p') = (q + p,p + i^sin(27rg' / )/(27r)), with kicking 
strengths K = 2.7 and K = 2.9, where one has a 



dominant partial barrier separating two sufficiently large 
chaotic regions, see Fig. [4] inset and Ref. p~4j [15] for de- 
tails. We consider the asymptotic transmitted weight Woo 
of a wave packet initially located outside of the partial 
barrier in region 1. We integrate the Husimi function of 
the wave packet ipi(t) over the measuring box £?2, which 
we choose as the entire region 2. Figure [4] shows good 
agreement of over the accessible range <&/h > - with 
the universal behavior as observed in Figs, [I] and [3] and 
well described by Eq. Q. 

We believe that the understanding of the universal be- 
havior of the quantum localizing transition of an isolated 
partial barrier provides the building block to explain the 
power-law scaling [5j [16] [17] occurring in the presence 
of hierarchically organized partial barriers around a can- 
torus. Our results also open the possibility to tackle the 
tunneling regime which sets in for <&/h < 0.1. There 
in addition to the local channel coupling mechanism of 
the turnstile one has tunneling across the entire partial 
barrier. Finally, if an extended chaotic system has an 
infinite chain of well isolated partial barriers, then the 
classical dynamics is diffusive, and the quantum dynam- 
ics will lead to exponential localization no matter how 
open the partial barriers. This similarity to Anderson lo- 
calization [9] would be very interesting to explore based 
on the universal transition curve of quantum transport 
through a single partial barrier. 

We are grateful to T. Guhr, M. Korber, J. Kuipers, 
and U. Kuhl for stimulating discussions and acknowl- 
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FIG. 4. (Color online) Asymptotic transmitted weight Woo 
using the Husimi measure for the standard map with K = 2.7 
($ « 250) and K = 2.9 ($ «_ ^) for ft = ^ . . . , ^ 
versus 4>/h compared to Eq. Q (solid line). The data is 
averaged over 15 initial conditions placed in region 1, 10 values 
of the Bloch phase, and 100 time steps after time 10 6 . Data 
for four different variants [141 115] of the partial barrier are 
shown. Inset: One partial barrier (thick solid line) separating 
regions 1 and 2 for the standard map with K = 2.7 as well as 
its preimage (dotted lines) 
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